Searching for the statistically equilibrated systems formed in heavy ion collisions 
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Further improvements and refinements are brought to the microcanonical multifragmentation 
model [Al. H. Raduta and Ad. R. Raduta, Phys. Rev. C 55, 1344 (1997); ibid. 61, 034611 (2000)]. 
The new version of the model is tested on the recently published experimental data concerning 
the Xe+Sn at 32 MeV/u and Gd+U at 36 MeV/u reactions. A remarkable good simultaneous 
reproduction of fragment size observables and kinematic observables is to be noticed. It is shown 
that the equilibrated source can be unambiguously identified. 
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I. INTRODUCTION 

It is known from more than 15 years that violent heavy 
ion collisions lead to an advanced disassembly of the com- 
pound system known under the name of nuclear frag- 
mentation. The asymptotically resulted fragments are 
the only experimental link with the primordial process. 
Good agreements between various observables related to 
these fragments and results of various models assuming 
uniform population of the system phase space |0, ^, |^ 
led to the conclusion that a statistically equilibrated nu- 
clear source is at the origin of the fragmentation process. 
The source size, its excitation energy and its volume are 
thus quantities which can only be indirectly evaluated 
by comparisons between experimental data and statis- 
tical multifragmentation models predictions. The com- 
parison process is complicated by various quantities pol- 
luting the pure statistical signals such as preequilibrium 
particle emission, collective radial expansion. Coulomb 
propagation of the break-up fragments, secondary parti- 
cle emissions. A complete theoretical reproduction of ex- 
perimental fragmentation event topology (never achieved 
so far) would help to answer questions fundamental for 
determining the thermodynamics taking place in such 
systems. For example, the identification of the experi- 
mental freeze-out volume would allow to locate the pro- 
cess in the phase diagram and, therefore, would help solv- 
ing a problem which made the subject of a strong debate 
over the last decade: is there any phase transition taking 
place in nuclear matter? 

The aim of the present paper is to investigate whether 
an unambiguous determination of the experimental sta- 
tistically equilibrated source can be achieved by means of 
the microcanonical model from Ref. |§, |^. To this aim, 
the model is refined by including the conservation of the 
angular momentum and completed with a Coulomb prop- 
agation stage of the primary fragments with the possibil- 
ity of superimposing radial flow. A complete set of ex- 
perimental fragmentation data concerning the reactions 
Xe+Sn at 32 MeV/u and Gd+V at 36 MeV/u recently 



published in Ref. Q is used for comparison with the 
model results. The paper is structured as follows: Sec- 
tion II gives a brief presentation of the employed micro- 
canonical model and describes its new refinements and 
improvements. In Section III a comparison between the 
model prediction and the above mentioned experimental 
data is presented. The influence of various model param- 
eters on different observables is discussed in Section IV. 
Conclusions are drawn in Section V. 



II. THE MODEL: BRIEF REVIEW AND 
IMPROVEMENTS 

The microcanonical multifragmentation model used in 
the present study was proposed in Ref. ||, |5| . The mod- 
ifications from the previous version of the model concern 
the primary break-up stage, momentum generation, in- 
clusion of collective radial expansion. Coulomb propaga- 
tion, and the secondary emission stage. 
i) Break-up stage 

This stage concerns the disassembly of a statistically 
equilibrated nuclear source (A, Z, E, V) (the mass num- 
ber, the atomic number, the excitation energy and the 
freeze-out volume respectively) . The basic assumption of 
the model is equal probability between all configurations 
C : {Ai, Zi, £i, Ti, Pi, i = 1,. . . ,N} (the mass num- 
ber, the atomic number, the excitation energy, the posi- 
tion and the momentum of each fragment i of the con- 
figuration C, composed of N fragments). Fragments are 
assumed to be spherical, are not allowed to overlap each 
other and are placed into a spherical recipient of volume 
V. In the previous version of the model |Q the system 
was considered subject to the standard microcanonical 
constraints: J^i = Yl,i '^-i = Z, J2i Pi = P (=0 in 
the cm. frame), E - constant. An extra constraint will 
be considered herein: the conservation of the angular mo- 
mentum, Ti X Pi = L. The integration over fragment 
momenta can be analytically performed subject to the 
above mentioned constraints: 
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Here Af = I is the inertial tensor of the sys- 

tem: Iaf3 = Y.i'^iij'i^ap - rfr'^^ with a,/? = 1,2,3, 

L^I-iL = i/ = E.P?/(2rn,:) + 

E.<, + E.e. - E.S. and = i? - E.<, - 
Ei + Ei V^ij stands for the Coulomb interaction 
between fragments i and j). Of course, the model cor- 
responds to the cm. frame (i.e. with i — 1,...,N 
are positions in the cm. frame) where P = [Q. If 
one also imposes L = (which is the hypothesis em- 
ployed in the present calculations) then the last factor 
m eq. (0) becomes (27r^C)3/2^-2. The only modifica- 
tions brought by the inclusion of the angular momentum 
conservation constraint to the statistical weight of a con- 
figuration C" : {Ai, Zi, Ei, Ti, z = 1, . . . , N} (W c) from 
Ref. 01 are the inclusion of the factor 1 / Vdet I and the 
replacement of N with — 1. These new weights can be 
employed in a Metropolis-type simulation which allows 
the determination of the average value of any system ob- 
servable in the very same manner as in Ref. B. 

Fragments with A < A are considered without excita- 
tion degrees of freedom except for the a particle for which 
the few levels larger than 20 MeV with T < 2.01 MeV 
have been considered. These fragments are weighted in 
Wc' with their energy levels degeneracies. Larger frag- 
ments carry internal excitation. For them, the following 
level density formula is included in the statistical weight 
of a configuration C (see Ref. [^): 

/tt 

= 12 ai/4e5/4 «^^P(2\/ae) exp(-e/r) (2) 

with a = 0.114A + 0.098^2/3 MeV^i § and r = 9 MeV. 
The factor exp(— e/r) is introduced to account for the 
dramatic decrease of the excited levels lifetime at high 
excitation energies 

ii) Generation of the primary decay fragment momenta 
While integration over the fragments' momenta has been 
carried out in order to simplify the Metropolis simulation, 
to produce events similar to the experimental ones one 
has to generate momenta for each given configuration C . 
Each of these events is characterized by a kinetic energy 
K. Therefore, we have to deal with the following com- 
putational task: generate uniformly momenta for a sys- 
tem composed of N particles such that J2iPi/'^''^i — 
Ei Pi = and X Pi = L. This problem was solved 

in an elegant way in Ref. Q . There it is proved that the 
following generation gives the right sampling of the sys- 
tem momenta, in agreement with the above mentioned 



constraints: Pick a preliminary set of N particle mo- 
menta from a arbitrary canonical distribution, then elim- 
inate the overall translational and rotational motion by 
making the transformation pi pi—mi{P'/M + uj' x r^) 
(where P' — p^ and a;' — J^i ^ Pi ' then spin 
the system such as to match the desired angular momen- 
tum L: Pi Pi + nii u) X Ti) with uj = h'^ ■ I^^ and 
finally, renormalize the momenta such as to match the 
available energy K — ^ u) ■ h. 
Hi) Radial flow 

After generating the fragment momenta corresponding 
to the primary decay, inclusion of nonequilibrium effects 
such as collective radial expansion (fiow) can be easily 
superimposed. Here we use the following parameteriza- 
tion for the flow velocity of fragment i: vj, = vo {ri/R)°', 
with vq = vo(ri/ri) and a a real number defining the fiow 
profile. After including the radial fiow the momentum of 
the «th fragment reads: pj ^ pj -I- to^v^. 

iv) Coulomb propagation 

After break-up, hot primary fragments are supposed 
to suffer an expansive propagation under their mutual 
Coulomb interaction. This stage can be easily simulated 
by integrating the corresponding Newtonian equations of 
motion. Integration has been carried out up to 500 fm/c, 
when Coulomb interaction between fragments can be ne- 
glected. 

v) Secondary decays 

Since primary fragments carry internal excitation, a sec- 
ondary decay stage was introduced in Ref. Q . Depend- 
ing on the fragment excitation, secondary break-up pro- 
cesses or particle evaporation were considered. Since the 
above classification is rather arbitrary, here we resume to 
the second process, treated using the standard Weisskopf 
evaporation scheme. As in the range of the evapo- 
rated particles is considered up to A = 16. Evaporation 
events are simulated using standard Monte Carlo [||. 

For each primary break-up event (i), the steps ii) - v) 
are performed. The resulted fragmentation events can be 
readily compared with experimental ones (after removing 
the particles coming from the preequilibrium stage). 



III. COMPARISON OF THE MODEL 
PREDICTIONS WITH RECENT 
EXPERIMENTAL DATA 

In this section comparison between the model results 
and the recent experimental data concerning the reac- 
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FIG. 1: Calculated charge multiplicities in comparison with 
experimental data for the Xe+Sn at 32 MeV/u and Gd+U at 
36 MeV/u reactions. Open circles correspond to calculation 
using the freeze-out hypothesis (1), open triangles correspond 
to calculation using the hypothesis (2) . Experimental data are 
represented by the histogram. 



FIG. 2: Calculated distributions of the largest, second 
largest and third largest charge in one fragmentation event 
in comparison with the corresponding experimental data for 
the Xe+Sn at 32 MeV/u and Gd+U at 36 MeV/u reactions. 
Symbols are used as in Fig. 1. 



tions Xe+Sn at 32 MeV/u and Gd+U at 36 MeV/u are 
presented. For the above reactions, various experimen- 
tal fragment size and kinematic observables have been 
pubhshed recently |^ . The aim of the study is to ex- 
tract information about the physical phenomenon taking 
place by trying to fit the entire set of experimental data 
using the microcanonical multifragmentation model. It 
is known [|lo| that various freeze-out assumptions may in- 
duce important differences in the statistical models' re- 
sults. For drastically reducing this uncertainty in the 
present study we employ two opposite freeze-out scenar- 
ios: 

1) The standard working hypothesis of the model: Frag- 
ments are idealized as hard spheres placed into a spherical 
freeze-out recipient; fragments are not allowed to overlap 
each other or the recipient wall. The generation of the 
fragments' positions is described in Ref. 

2) The hardcore interaction is switched off: Integration 
over the fragments' positions may be approximately car- 
ried out by assuming as in Ref. that each fragment 
(from a configuration composed of N fragments) is block- 
ing the volume Vq/N {Vq is the volume of the nuclear 
system at normal density) for the rest of the fragments 
(as in [|l2j) and for itself as well. Coulomb energy is ap- 
proximated by the Wigner Seitz formula, being therefore 
independent of the fragments' positions Oq]. The inte- 



gration over the fragments' positions writes: 
/■ ^ 1 r ^ 1 
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The factor Y\j^i Oij from eq. (||) is just a formal expres- 
sion of the fragment blocking constraint formulated be- 
fore. With this approximation fragment positions have to 
be generated into the spherical freeze-out volume, with- 
out any hardcore constraint only for evaluating det I and 
the statistical weight of a configuration will get (as in 
Ref. ||l^) the extra factor Vfree- Having the weights of 
each configuration C" the corresponding Metropolis type 
simulation is straight forward |^|. 

Using the above defined working hypotheses, compar- 
isons between the model's results and the corresponding 
experimental data have been performed. Both fragment 
size distributions and kinematic distributions have been 
considered. Very good agreement between calculations 
and experimental data can be observed for both reac- 
tions, in both considered working hypotheses for all the 
considered observables. 
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FIG. 3: Calculated IMF multiplicities (left column) and 
Zbound3 multiplicities (right column) in comparison with the 
corresponding experimental data for the reactions Xe+Sn at 
32 MeV/u (upper raw) and Gd+U at 36 MeV/u (lower raw). 
Symbols are as in Fig. 1. 

The comparisons between theoretical and experimental 
fragment size distributions corresponding to both Xe+Sn 
at 32 MeV/u and Gd+U at 36 MeV/u reactions are as 
follows: Charge multiplicities are given in Fig. 1. Multi- 
plicities of largest charge (Zmaxi), second largest charge 
{Zmax2) and third largest charge (Zmaxs) from one frag- 
mentation event are presented in Fig. 2. In Fig. 3 inter- 
mediate mass fragment (IMF) distributions and Zbounds 
distributions are given (IMF fragments are considered to 
have Z > 5 as in 01 ; Zbounds is defined as the sum over all 
fragments' atomic numbers from one fragmentation event 
which are greater or equal to 3; Zhoundb is similarly de- 
fined but the Z limit is fixed to 5) . Zhoundb multiplicities 
are compared with experimental data in Fig. 4 only for 
the reaction Gd+U at 36 MeV/u. It can be noticed that 
for the Zbound distributions the calculated width of the 
distributions are slightly smaller than the experimental 
ones. This suggests a (small) fluctuation of the experi- 
mental equilibrated source. This will become more clear 
in the next section. An excellent agreement for all consid- 
ered observables related to fragment size was, therefore, 
obtained. 

Comparisons between theoretical and experimental av- 
erage fragment kinetic energy function of fragment charge 
are given in Fig. 5, left column for the the two consid- 
ered reactions. Very good agreements are to be noticed. 
A deeper insight to the kinematics of the process can 
be obtained by analyzing the reduced fragments velocity 
correlations ||ll[. Fragments with 5 < Z < 20 as in fl^ 
have been used for constructing the correlations. This 
observable was proven to reflect subtle topological fea- 
tures of the fragmentation events. The beginning of the 



correlation function is often called "Coulomb hole" and 
reflects the strength of the Coulomb repulsion between 
fragments (obviously related to the size of the freeze- 
out volume) but, as shown in the next section, is also 
related to the amount of radial expansion. Good agree- 
ment between calculations and experimental data can be 
observed here as well (see Fig. 5). While in the case 
of the Gd+U reaction the experimental velocity corre- 
lations are well reproduced in both freeze-out hypothe- 
ses, in the case of the flrst reaction a better agreement 
with experimental data is obtained for the hypothesis 
(2) ("without hard-core"). In fact, this result is natural 
and suggests that deviations from the spherical fragments 
with hard-core interaction idealization (deformed frag- 
ments, surface diffusivity, etc.) are present in real mul- 
tifragmentation. The deviation effect is stronger in the 
case of the Xe+Sn since the corresponding equilibrated 
system is smaller and, therefore the fragment partitions 
are more symmetric (as shown in the free volume 
is smaller for more symmetric partitions and larger for 
more asymmetric ones). The simultaneous description of 
both average fragment kinetic energies as a function of 
the fragment charge and of the reduced velocity correla- 
tions is remarkable and reflects a good description of the 
event topology. Here it is worth noticing that in Ref. Q 
the velocity correlations for both reactions could not be 
reproduced via the SMM model. 

The values of the model input parameters for which 
the comparisons with the experimental data have been 
performed are listed in Table |. A remarkable fact is 
that the fltting parameters corresponding to the freeze- 
out hypotheses (1) and (2) are quite close ones to the 
others for each of the considered reactions. This tends 
to indicate a relatively small "region", speciflc to each 
reaction in which the parameters of the real equilibrated 
nuclear source are situated. Quite small differences from 
one hypothesis to the other seem to be present in the case 
of the freeze-out volume. This result may look intriguing 
given the fact that freeze-out volume is directly related 
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FIG. 4: Same as in Fig. 3 (right column) for the Zbounds 
distribution corresponding to the Gd+U at 36 MeV/u 
reaction. 
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FIG. 5: Average fragment kinetic energy versus charge (left 
column) and reduced fragment velocity correlation (right 
column) for the reactions Xe+Sn at 32 MeV/u (upper raw) 
and Gd+U at 36 MeV/u (lower raw). Symbols are as in 
Fig. 1 except for the experimental data which are here 
represented by full circles. 

to both working hypotheses. In reality, the relatively 
large freeze-out volumes obtained are minimizing the in- 
fluence of the hard-core interactions to the free-volume 
(Yfree) (sec e.g. ) , the values of of the freeze-out vol- 
ume being mainly dictated by the Coulomb interaction 
(the remaining quantity entering the system's density of 
states depending on the freeze-out volume). Unlike the 
free volume which in models can only be more or less 
arbitrary parametrized, Coulomb interaction can be ac- 
curately evaluated for each fragment partition (either by 




TABLE I: Values of the model's input parameters corre- 
sponding to the fitting of the experimental data for both 
Xe+Sn at 32 MeV/u and Gd+U at 36 MeV/u reactions 
in the freeze-out hypotheses (1) and (2). The last two 
columns correspond respectively to the flow energy and to 
the exponent defining the flow profile (see text). 



react. /hyp. 


{A,Z) 


V/Vo 


£;,4MeV/u) 


£;/,(MeV/u) 


a 


Xe+Sn (1) 


(230,94) 


9 


5.3 


1.4 


2 


Xe+Sn (2) 


(220,92) 


8.5 


4.78 


1.9 


1.8 


Gd+U (1) 


(343,136) 


8.5 


5.95 


1.7 


2 


Gd+U (2) 


(328,130) 


8 


5.2 


1.9 


1.8 



FIG. 6: Zboundi distributions (symbols) corresponding 
to various sizes of the nuclear source within the working 
hypothesis (1) (see the legend). 



explicitly accounting for the interfragment interaction, 
or by employing the Wigner-Seitz approach). This con- 
tributes to the solidity of the present evaluations. 

The differences in the excitation energies correspond- 
ing to the two hypotheses are due to the Coulomb en- 
ergy variation from one case to another, caused by the 
slightly different values of the freeze-out volumes (see Ta- 
bic H). For example, in the Gd+U case the energy varia- 
tion (from (1) to (2)) is slightly larger than in the Xe+Sn 
case since the relative variation of the freeze out volume 
is larger as well. Freeze-out volume is always smaller in 
the working hypothesis (2). This is related to the differ- 
ent dependencies VfreeiY) resulting from the underlying 
hypotheses (more symmetric partitions and thus a more 
advanced fragmentation are specific to the without hard- 
core case and therefore both the freeze-out volume and 
excitation energies need to be smaller in the second hy- 
pothesis). The flow parameters were chosen such as to 
insure a good reproduction of the kinetic observables (see 
Fig. 5). The way in which each parameter is influencing 
the fitting observables is discussed in the next section. 



IV. INFLUENCE OF THE MODEL 
PARAMETERS ON VARIOUS FITTING 
OBSERVABLES 

While many fragment size and kinetic observables were 
simultaneously fitted in order to deduce the parameters 
of the equilibrated sources corresponding to the two re- 
actions, the following question arises: are they unique? 
We address this question in the present section. 

We start with the remark that the size of the equili- 
brated source appears to be dictated by the Ztound dis- 
tribution. We illustrate this by representing in Fig. 6 
the Zhound3 distributions corresponding to three differ- 
ent sources: (215,90), (230,94) and (245,102). Their ex- 
citation energies (5.4, 5.3, 5.3 MeV/u respectively) are 
chosen such as to insure a good fit of the charge distribu- 
tions corresponding to the Xe+Sn at 32 MeV/u reaction; 
their freeze-out volume {V/Vq = 9) is chosen as to pro- 
vide a good description of the Z^ax, Zmax2 and Zmax3 
distributions. It can be easily observed that the best fit 
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FIG. 7: Largest, second largest and third largest charge distributions corresponding to two different freeze-out volumes (see 
the legend) for each of the two freeze-out hypotheses ((1) upper raw; (2) lower raw) in comparison with Xe-|-Sn at 32 MeV/u 
the experimental data. See the text for the sources' parameters. 



on the data is given by the (230,94) source. A smaller 
source deviates the maximum of the Zhounds distribution 
to the left while a larger one deviates it to the right (see 
Fig. 6). 

The freeze-out volume of the equilibrated sources ap- 
pears to be dictated by the distributions of the largest, 
second largest and third largest charge in one fragmenta- 
tion event. Obviously, these distributions provide a mea- 
sure of the degree of asymmetry of a given fragment parti- 
tion. The degree of asymmetry of a partition is influenced 
by the freeze-out volume through both the Coulomb in- 
teraction and the free volume . This fact is evidenced 
in Fig. 7 where multiplicities of Zmaxi, Zmax2 and Z^axz 
are represented for two different freeze out volumes in the 
freeze-out hypotheses (1) and (2) as follows. For the first 
hypothesis the source was considered (230,94) at two dif- 
ferent volumes: Vb/5.5 and Vb/9 while for the second one 
the source was taken to be (220,92) at other two different 
volumes Vb/3 and Vo/8.5. The corresponding excitation 
energies are taken such as a good reproduction of the 
experimental charge distributions corresponding to the 
Xe-|-Sn reaction at 32 MeV/u to be achieved (the excita- 
tions energies are respectively: 6.4, 5.3, 6.1, 4.78 MeV/u). 
It can be noticed that deviations to the left of the max- 
ima (indicating more asymmetric fragment partitions) of 
the Zraaxi and Zmax?, distributions are occurring when 
the freeze-out volume is decreased. This behavior can 
be observed irrespective to the considered freeze-out hy- 
pothesis or to the chosen mass of the source. While, as 
demonstrated by Fig. 7 lower values of the freeze-out vol- 
ume are leading to deviations of the Zmaxi, Z^aaxi and 
Zmax?, calculated distributions from the corresponding 
experimental data, larger values of the freeze-out volume 



will result in smaller excitation energies and larger flow 
energies leading to deviations from the experimental data 
of the kinematic observables represented in Fig. 5. 

Given the previously described monotonical behavior 
of the "distance" from the calculated curves to the ex- 
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FIG. 8: Influence of the radial flow parameters on the frag- 
ment average kinetic energy versus charge dependency and 
on the fragment reduced velocity correlations for the working 
hypothesis (1) (see the legend and the text for the sources' 
parameters). The experimental data are here represented by 
full circles. 
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perimental ones, and also the very good reproduction of 
all the considered experimental observables for the pa- 
rameters (given in Table I) corresponding to the case in 
which the Ztounds distribution is fitted best (see Fig. 6) 
it results that the corresponding sets of parameters pro- 
vide the global minimum of the associated function. 
Indeed, choosing a different value of A would generate a 
deviation of the Zbounds from the experimental one and, 
obviously the global reproduction of the data would be 
worst. 

Finally, the way in which the radial flow parameters 
Efi and a are influencing the experimental data fitting 
is presented in Fig. 8 for the nuclear source (230,94), 
V = Vo/9, Eex = 5.3 MeV/u and the first working freeze- 
out hypothesis. Three different cases are being consid- 
ered: [a = 2, Efi = 1.4 MeV/u), (a = 1, Efi = 1.4 
MeV/u) and [a = 2, Eji = 0.7 MeV/u). As shown in 
Fig. 8 the a parameter defining the flow profile appears 
to influence the location of the maximum of the aver- 
age kinetic energy versus Z {Ec.m.{Z)) and the value 
of the flow energy the height of the Ec.m.i^) depen- 
dency. Smaller effects are to be noticed in the reduced 
velocity correlation function and they concern mainly 
the width of the "Coulomb hole". The effects are in- 
tuitive: smaller flow leads to a narrower Coulomb hole 
while larger ones lead to a larger hole. The parameters 
a = 2 and Efi — 1.4 MeV/u appear to give the best de- 
scription of the kinetic experimental data corresponding 
to the Xe-|-Sn at 36 MeV/u reaction. 

Therefore, it results that, within some (small) inherent 
uncertainties, the obtained parameters corresponding to 
the two working hypotheses ((1) and (2)) are uniquely 
determined. 



V. CONCLUSIONS 

In summary, the microcanonical multifragmentation 
model 1^ has been further refined and improved. Con- 
servation of angular momentum has been included. Also, 
in order to provide events similar with the experimental 
ones, the model is completed with a stage of generation 
of momenta for the primary break-up fragments. Radial 
flow velocities can be superimposed here. Finally, an 
expansion stage of the hot primary fragments under 
their mutual Coulomb interaction is added. Using the 
new version of the model, the recent experimental data 
of the INDRA collaboration concerning the reactions 



Xe+Sn at 32 MeV/u and Gd+\J at 36 MeV/u are 
analyzed. Both fragment size distribution observables 
and kinematic observables could be simultaneously 
fitted allowing thus the identification of the equilibrated 
sources parameters corresponding to the two reactions. 
Two different break-up scenarios were used for perform- 
ing the analysis. The two types of calculations give close 
results for each of the considered reactions suggesting 
independence of the obtained results on the particular 
hypothesis. This is due to the fact that for relatively 
large values of the freeze-out volume (such as those 
obtained by fitting the data, 8-9 Vq) the influence of 
the hard-core is small and thus the results are rather 
independent on the extent to which these effects are 
included in a specific Vfree parametrization. In this 
case, the values of the freeze-out volume are dictated 
by the Coulomb interaction (accurately evaluated 
independently to any model assumptions) which makes 
the evaluation particularly robust. For each reaction a 
narrow freeze-out volume region was identified (8.5-9 Vq 
for Xe-hSn and 8-8.5 Vq for Gd-hU). By analyzing the 
way in which variations of the model parameters are 
influencing the deviation between the calculated curves 
and the experimental ones, it is proved that the set of the 
parameters providing the best fit to the data is unique 
(within inherent uncertainties). It is shown that the 
size of the source is dictated by the Zbound distributions 
while its freeze-out radius can be deduced by fitting 
the experimental distributions of Zmaxi, Z^axi and 
Zmaxz and the kinematic observables. The way in which 
the fiow energy and the parameter a, defining the fiow 
profile is influencing the fragment average kinetic energy 
versus Z and also the reduced velocity correlation is 
also discussed. The above parameters are shown to be 
clearly identifiable as well. To our knowledge, this is the 
first complete reproduction of an entire set of experi- 
mental data (including both fragment size distribution 
observables and kinematic observables) by means of a 
microcanonical multifragmentation model. Moreover, 
it was shown for the first time that the parameters of 
the experimentally obtained equilibrated source can be 
uniquely determined by means of such a model. 
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